WIMP Dark Matter Inflation with Observable Gravity Waves 
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We present a successful inflation model based on Xcj)'^ potential in which a Standard Model (SM) singlet 
inflaton </!>, with mass of around a TeV or less, also plays the role of a weakly interacting scalar dark matter 
particle (WIMP). The WIMP relic abundance generated after inflation is in accord with the current observations. 
The spectral index ris lies within the WMAP I-ct bounds, while the Planck satellite may observe the tensor- 
to-scalar ratio, a canonical measure of gravity waves, which we estimate lies between 0.003 and 0.007. An 
unbroken Z2 parity ensures that the scalar WIMP is absolutely stable. 
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The idea that the inflaton, a particle responsible for primor- 
dial inflation, also may play the role of scalar WIMP dark 
matter is most intriguing |1] and therefore worth pursuing. 
Ref. ^ attempted to implement this idea in chaotic inflation 
with m^(/)^ potential. However, a satisfactory scenario could 
not be realized which is, to a large extent, related to the fact 
that m ~ 10^'^ GeV, as demanded by inflation, far exceeds the 
canonical WIMP mass of a TeV or so. Simply replacing the 
quadratic potential with a quartic one does not help solve the 
conundrum for in this case the scalar spectral index Ug and 
tensor-to-scalar ratio r lie outside the WMAP 2-a bounds ||3|]. 

In a recent paper, hereafter called [4], it was shown that 
A(/)^ inflation, if supplemented by the non-minimal gravita- 
tional coupling ^TZc/P between the SM gauge singlet scalar 
field <j) and the curvature scalar TZ, yields values of Ug (scalar 
spectral index) and r (tensor-to-scalar ratio) that are compati- 
ble with the WMAP 1-cr bounds yl]. This is to be contrasted 
with minimal Xcf)'^ inflation whose predictions for n^ and ?■ 
lie outside the WMAP 2-a bounds. A series of earlier papers 
OSll- illll have previously raised the possibility that the inflaton 
field (p could be identified with the SM Higgs doublet H, pro- 
vided the non-minimal coupling ^ = C9(10^)-C9(10*). While 
intriguing, doubt about the viability of this identification have 
been raised in 1 1211 il3ll . It stems from the observation that for 
^ ^ 1, the energy scale X^^'^nip/y^ of inflation exceeds the 
effective ultraviolet cutoff scale mp/£^, with mp being the re- 
duced Planck scale, assuming the SM Higgs quartic coupling 
A is of order unity. In [4] we easily evade this problem by 
making (j) a SM gauge singlet field so that the parameter A is 
not all that strongly constrained. Indeed, one finds that con- 
sistent with the WMAP 1-cr bounds on n^ and r, A and ^ can 
lie within the relatively wide range, 10^^^ ^ ^ S 10^'', and 

In a separate development, it has been noted by several au- 
thors [ 14] that a stable SM singlet scalar particle, with mass ^ 
m/i/2— 1 TeV, is a viable cold dark matter candidate (WIMP), 
provided it has suitable interactions with the SM Higgs dou- 
blet H and possibly additional fields. The interaction term 
g^(fP\H\^ plays an especially important role in these consid- 
erations. Recent estimates suggest Iil5i] that with g^ ~ 0.1 
and dark matter mass ^ 1 TeV, the relic WIMP abundance is 
compatible with the value ricDM/i^ = 0.1131±0.0034deter- 
mined by WMAP liSJ. This parameter region will be further 



explored in the ongoing and planned direct detection experi- 
ments of dark matter particle. 

In this letter we propose a successful and relatively sim- 
ple scenario of WIMP dark matter inflation by merging to- 
gether ideas from ||4J] and lIlSll . Following |4], we employ non- 
minimal quartic inflation in which a gravitational coupling of 
the inflaton to the curvature scalar is included. The model has 
a further restriction arising from the relic dark matter abun- 
dance. It is shown in [ISil that for TeV mass WIMP dark mat- 
ter, the coupling strength g'^ must be of order 0. 1 or so. In 
our case this means that due to radiative corrections involv- 
ing g^, the 'effective' quartic coupling is of order 10^'^. An 
important consequence of this WIMP driven inflation model 
is that it predicts both Ug and r in a fairly narrow range. In 
particular, r values close to 0.007 may be accessible to Planck 
satellite searches. Another important feature of our model is 
the appearance of thermal dark matter relic abundance which 
arises during preheating and subsequent transition to a radia- 
tion dominated universe with temperature close to 10^ GeV. 
The energy in the oscillating inflaton field is by then essen- 
tially negligible. 

Consider the following tree level action in the Jordan frame: 
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where rnp ~ 2.4 x 10^^ GeV is the reduced Planck mass. 
Here we have introduced a Z2 parity under which (j> is odd, 
while the SM fields are all even. Hence the scalar (p is stable 
and will play the role of both inflaton and dark matter particle. 
First we consider the non-minimal X(j>* inflation in this 
model |4]. During inflation, with field values close to mp, 
A0* dominates the scalar potential. The relevant one-loop 
renormalization group improved effective potential 11611 is 
Veff = iA(i)G(i)V'^, where t = ln((/i/m0), and G{t) = 
exp(- J*dt'j{t')/{1 + 7(i')))' with j{t) being the anoma- 
lous dimension of the inflaton field. We employ a leading-log 
approximation for the effective potential 
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where Aq = X{t = 0), Aq ^ 5^, and we have taken m^ as the 
renormahzation scale. In the Einstein frame with a canonical 
gravity sector, the kinetic energy of (p can be made canonical 
by defining a new field a I?]], 



da\ 



1 



?</>' 



(3) 



The effective potential in the Einstein frame is then given by 



The inflationary slow-roll parameters are given by 
1 . 



(4) 
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where a prime denotes a derivative with respect to </>. The 
slow-roll approximation is valid as long as the conditions e <C 
1, I77I ^ 1 and C,^ <^ 1 hold. In this case the scalar spectral 
index Ug, the tensor-to-scalar ratio r, and the running of the 
spectral index a — jj^ are approximately given by 
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The number of e-folds after the comoving scale I has crossed 
the horizon is given by 
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where 0/ is the field value at the comoving scale I, and 0e 
denotes the value of cj) at the end of inflation, defined by 
max(e((/)e), |?7((?!)e)|, C^('/'e)) = 1- The amplitude of the cur- 
vature perturbation A-ji is given by 
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which should satisfy the WMAP normalization, A^ = 
(2.43 ± 0.11) X 10-^ |3], at fco = 0.002 Mpc"\ 

Using Eqs. (|2]l-(l8]l we can obtain various predictions of the 
radiatively corrected non-minimal XcfA inflation model. Once 
we fix the parameters ^ and the number of e-foldings A^q, we 
can predict n^, r, and a — jj^- Note that with m^ ~ 1 
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FIG. 1: r vs. Us with No — 60 (solid curve) and A'^o — 50 (dashed 
curve) e-foldings. Both curves lie within the WMAP 1-a (68% con- 
fidence level) bounds. 



TeV, the coupling g^ ~ 0.1 in order for the relic density of 
dark matter to be compatible with the WMAP observations 
OlSll . In our analysis, we set gi^ = 0.1 and m^ = 1 TeV 
as reference values. If Aq < 5^/(87r'^), the potential during 
inflation is dominated by the radiatively corrected part. We 
impose Aq > for an unbroken Z2 parity. 

The predicted values of Ug and r are shown in Figure [T] 
for the number of e-foldings A'o — 60 (solid curve) and 
A^o — 50 (dashed curve). In non-minimal Xcj)'^ inflation ||4|], 
the (effective) scalar quartic coupling becomes larger accord- 
ing to ^ values, and Ug and r approach their asymptotic values. 
Us ~ 0.968 and r ~ 0.0030 for A^o = 60 and n^ ~ 0.962 and 
r ~ 0.0042 for Aq = 50. These values correspond to the 
left edge of each curve in Figure [1] In the present case, for 



Ao < g /(Btt ), the radiatively induced term in the effective 
potential dominates the scalar potential and thus the effective 
quartic coupling has a minimum value. In the limit Aq — 0, 
Us and r approach rig ~ 0.976 and r ~ 0.0054 for A^o — 60 
(n, ~ 0.970 and r ~ 0.0069 for A^o = 50), which correspond 
to the right edge of each curve in Figure [T] We find that the 
running of the spectral index a — jf^ very weakly depends 
on Us, and a ~ -0.0005 (-0.00075)for A^q = 60 (A^o = 50). 
In Figure |2] we show the ratio of the inflation energy scale 
(\/i/4) to the effective ultraviolet cutoff scale (A = rnp/S). 
This ratio becomes larger as ^ is raised. The minimum value 
for this ratio, V^/*/A ~ 8.7, is achieved for Ao < 5'*/(87r2) 
with ^ ~ 2000. This value marginally exceeds the proposed 
naturalness bound V^^'^/A < 0{1) 1 1 311 . As a more con- 
servative bound, we examine the constraint on the ratio of 
the Hubble parameter and the effective cutoff scale from the 
validity of the classical inflationary treatments 01211 . namely, 
VA < H/A < 1. Since H/A ~ {V^/'^/A)'^{A/mp), we 
find the ratio - (8 .7)V2000 ~ 0.04 for ^ ^ 2000, while 
VA ~ ^yg^/{8^^^) ~ 0.01. Thus, this bound is satisfied. 

After inflation, the inflaton field starts to oscillate around its 
potential minimum. During this period the energy density of 
the inflaton is transmitted to other relativistic fields. Since a 
single inflaton cannot decay into other particles because of Z2 
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FIG. 2: y^/^/A vs. n^ with No = 60 (solid curve) and A^o = 50 
(dashed curve) e-foldings. 



parity, preheating yj] via, in our case, the coupling g'^(jP\H\'^ 
plays the crucial role in energy transmission. In the first stage 
of preheating, the inflaton energy density is transmitted to 
(j> particles by explosive production through parametric res- 
onance effects with A0^ potential. The amplitude of the in- 
flaton filed is reducing its amplitude by particle production 
and the expansion of the universe. When this becomes smaller 
than TTi^/vA, the term ra^cjp' in the scalar potential domi- 
nates and the Higgs doublets are explosively produced by the 
(broad) parametric resonance through the coupling g^ (f)'^\H\'^ . 
This preheating process ends when the amplitude becomes 
smaller than in^/g lV\. 

Since the inflaton is stable, its oscillations continue without 
further energy transfer after the preheating era. The oscillating 
mode has the equation of state of dust and so, in principle, 
it can play the role of dark matter. The present ratio of the 
average energy density of this oscillation mode to the number 
density of photons is estimated as |2] 
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where 4>* is the oscillation amplitude at time i, when to^ = 
iJ*, the Hubble parameter Since m^ > H for t > t^,, (p^ 
is basically also the amplitude of the oscillating inflaton to- 
day. With 0, = m,cf,/g, we find ^dm,o — 1-5 x 10~^^top for 
m^ = 1 TeV and g^ = 0.1. Comparing it to the observed 
value by WMAP [3], ^dm,o — 1-1 x lO^^^mp, we conclude 
that the oscillating mode has a negligible contribution to the 
energy density of the present universe. It is worth recalling 
that successful m|(/)^ chaotic inflaton requires m^ ~ lO^'^ 
GeV, and so g ~ lO'^ is needed to realize the observed value 
of '?dm,o B2ll! In our case, successful inflation is realized by 
non-minimal X(j>'^ inflation and the inflaton mass plays no role 
during inflation for m^ ^ lO^'^ GeV. 

We have seen that the energy density in the remnant inflaton 
oscillations is tiny compared to the observed dark matter relic 
density. However, the 4> particle can be a suitable WIMP dark 
matter candidate if the universe is thermalized with the reheat- 



ing temperature high enough for the cj) particle to be in ther- 
mal equilibrium. In the preheating scenario, thermalization of 
the universe takes place through decays and multiple scatter- 
ings of particles (SM Higgs doublets in our model), during 
explosively produced preheating. A reasonable estimate for 
the reheating temperature is 1 1 ] 
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where Th is the total decay width of the Higgs boson. For a 
relatively light Higgs boson with mass ruh = 120 GeV for 
example, the dominant decay mode is /i — > bb, so that 
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where mt — 3 GeV is the appropriate bottom quark mass, and 
V = 246 GeV is the VEV of the SM Higgs doublet. We find 
Tp ^ 10^ GeV, which is four orders of magnitude larger than 
TO^(~ 1 TeV). Thus, we expect that the temperature of the 
universe is high enough for (p particles to be in thermal equi- 
librium, and the standard WIMP dark matter scenario consis- 
tent with the WMAP observations can be realized, as shown 
in recent analysis [15|]. 

In our analysis above, we set m^ = 1 TeV in order to keep 
the inflaton mass much larger than the SM Higgs boson mass, 
7710 ^ mil. This parameter choice makes the preheating pro- 
cess as effective as possible. In general, the preheating pro- 
cess can be reasonably efficient even for m^ ^ m.h 111]. If 
the inflaton mass can be lowered close to half of the Higgs 
boson mass, the magnitude of the coupling g^ needed to ob- 
tain the correct dark matter relic abundance becomes signif- 
icantly smaller than 0.1 PIS']. In this case, the radiative cor- 
rections to the potential is negligible and we can easily obtain 
V^l^ j tv < 1 for a successful inflation scenario as shown in 
li4J. 

In summary, we have shown that the inflaton and WIMP 
dark matter can indeed be one and the same particle. In the 
simplest model this is achieved by supplementing the SM with 
a stable gauge singlet scalar field. The model overcomes se- 
rious challenges faced by chaotic 777^0^ inflation |2] and, in 
addition, turns out be quite predictive. Its dark matter prop- 
erties will be seriously examined by the ongoing direct detec- 
tion searches lllTll lllSll . As far as inflation is concerned the 
predictions for 77s and r lie within the WMAP l-a bounds. 
With an upper bound of around 0.007 on r, the model can be 
excluded if the Planck satellite observes values that are signif- 
icantly larger than this. 

Finally, one promising extension of our model is to intro- 
duce an SU(2) triplet scalar field with unit SM hypercharge. 
This would nicely incorporate neutrino masses and mixings 
via the type-II seesaw mechanism il9ll . Interestingly, this ex- 
tension essentially coincides with a model proposed in |2C|], 
and it can also account for the anomalous cosmic -ray positron 
flux reported by the PAMELA satellite experiment 112 ill . In 
addition, as analyzed in detail in II22I1 . in the type-II seesaw 
extension of the SM, the vacuum stability bound on the Higgs 



boson mass can be reduced to coincide with the current exper- 
imental lower bound of 1 14.4 GeV ll23ll . 



Note Added 

Although our approach in unifying the inflaton and dark mat- 
ter particle is inspired by Refs. yj] l|2|l, by the recent analysis 
of SM singlet^calar dark matter lllSn . and by non-minimal 
A(/>^ inflation f4], the model presented in this letter turns out 
to be identical to the ones proposed in Refs. 1241 and [25]. In 
II24II . a comprehensive study of certain aspects of this model 
were presented, and where there is overlap with our work, the 
results appear to be in broad agreement. However, we em- 
phasize that there are several new and important results in this 
letter Thus, (1) we have emphasized the WMAP constraints 
on the coupling g^, arising from the thermal relic density of 
dark matter, and evaluated its direct impact on the effective in- 
flaton quartic coupling. As a result, our inflationary scenario 
predicts both n^ and r in a narrow range. (2) This feature also 
then plays an important role in discussing the naturalness of 
the inflationary scenario as shown in Figure |2] (3) We have 
considered the preheating scenario after inflation, following 
Ul and II2I], which allows the transition to a radiation domi- 
nated universe. Furthermore, we have shown that the energy 
density in the remnant inflaton oscillations, which plays the 
role of dark matter in the original scenario ||2|] , can be ignored. 
(4) Considering thermalization of the universe via preheating, 
we have estimated the reheating temperature and shown that 
it is high enough for the singlet scalar to be in thermal equi- 
librium. This is crucial for a successful WIMP dark matter 
scenario. 
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